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Abstract 

Topologically nontrivial time-dependent solutions of the classical non- 
linear a model are studied as candidates of the disoriented chiral con- 
densate (DCC) in 3+1 dimensions. Unlike the analytic solutions so far 
discussed, these solutions cannot be transformed into isospin-uniform ones 
by chiral rotations. If they are produced as DCCs, they can be detected 
by a distinct pattern in the angle-rapidity distribution of the neutral-to- 
charged pion ratio. 
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1 Introduction 



Analytic solutions of the classical nonlinear a model [p], @, || have been studied 
as candidates of the disoriented chiral condensates (DCCs) ||. The solutions so 
far obtained are either configurations with spatially uniform isospin distribution 
or those which are chirally equivalent to them. When the isospin-uniform DCCs 
decay, the decay pions will obey the event-by-event pion charge distribution 
of dP/df = l/2v7 m the neutral pion fraction /. In this paper we study as 
DCC candidates the time-dependent solutions of the nonlinear cr-model that are 
topologically nontrivial in the isospin-orbital space. Though the topologically 
nontrivial DCCs obey the same charge distribution dP/df as that of statistically 
random emission, the angle-rapidity distribution of pions should exhibit a very 
distinct experimental signature. We suggest a quantitative method of analysis 
to search for these DCCs. 



2 Topologically nontrivial solutions 

Let us express the pion field ir{x) nonlinearly in terms of the scalar phase field 
9{x) and the unit isovector field n(x) as 

7r(x) = f„n(x)9(x), (1) 

where f n = 93 MeV is the pion decay constant. Apart from the source term, 
the Lagrangian is given in terms of = exp(m(x) -r 9(x)) by 

£ = S-tr(d^d^) + \\f^-l), (2) 

in the chiral symmetry limit. After elimination of the Lagrange multiplier A, 
the Euler-Lagrange equation reads @ 

□0 - sin 6 cos 6 (<9 M n • d^n) = 0, (3) 

<9 M (sin 2 # n x ^n) = 0. (4) 

The analytic solutions so far known are either those with n(x) = constant 
(the Anselm class) or those which are rotated to them by chiral transformations. 
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In this paper we explore the class of solutions whose n(x) fields point radially 
in the spacetime of 3 + 1 dimensions: 

n(x) = *-. (5) 

We have obtained a hint for this postulate from the Skyrme assumption that 
led to the static soliton |J. Unlike Skyrme, we do not need a stabilizing term 
in the Lagrangian since the static stability of solutions is irrelevant to us. The 
spherical symmetry of n(x) suggests we should choose 6{x) also to be spherically 
symmetric as 

e(x) = e(t,r). (6) 

It is important to observe that with our postulate the Euler-Lagrange equation 
(4) for n(x) is automatically satisfied for any 8(t, r). The equation for 6{x) now 
reads 

{d 2 t - V 2 )9(t, r) + \ sin 6(t, r) cos 6(t, r) = 0. (7) 

This wave equation allows many interesting solutions. All of them are topologi- 
cally nontrivial since following the Skyrme model we can introduce the topolog- 
ical charge, 

q = J Q d 3 x, (8) 

where 

= ^tr(X'X^), (9) 

with X u = E^c^E. The current is locally conserved, d^Q^ = 0, and q is 
invariant under chiral SU(2) x SU(2) rotations. The charge q is nonvanishing 
for our solutions while it is zero for the isospin uniform solutions. It should be 
noted that q is time dependent when we compute it for the pion fields alone since 
the current flows from the shell of hadron debris into the DCC Actually it 
is not even finite since our pion fields are singular as we approach the light cone. 
This should not bother us since the nonlinear a model after all does not apply 
to the close neighborhood of the light cone where kinetic energy is too large. 
Our purpose of mentioning the topological charge q here is that our solutions 
are chirally inequivalent to the isospin-uniform solutions. When we adopt our 
solutions as DCC candidates, we do not accept the Skyrme model of baryons 
0] in which the charge q is identified with the baryon number. If we did, our 
DCCs would be loaded with nucleons or antinucleons. 
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Let us solve for 9(t,r) by restricting the form of 9(t,r). Since 9 is dimen- 
sionless and its equation of motion is scale invariant, a simple case of interest is 
that 9(t, r) is a function only of the ratio of r and t: 

9(t,r) = 9(0, (10) 
The wave equation then turns into 

l^ + ^-^^sintfcos^O. (11) 

It can be easily solved numerically. The behavior at £ = is determined by the 



singularity at £ = in the wave equation (]TT|) . Barring a singularity at r = 
for 9(x) since a source does not exists at r = after t = 0, we determine that 
9(0 oc £ as £ — > 0. By giving one more boundary condition, d9jd^ at £ = 0, we 
can compute a profile of the scalar phase field In Fig.l we have plotted 9(£) 
for a few different values of 9 f (0). On the light cone (£ = 1), 9(t,r) is singular 
because of 1/(1 — £ 2 ) in the third term of Eq.flll]), but only in the immediate 
neighborhood of £ = 1. The function is smooth and monotonic practically 
everywhere inside the light cone. 

When 9(t, r) is not a function only of £ = r/t, analyticity at r = requires 
the behavior 9(t, r) — > r at r — > 0, not necessarily #(£, r) — > r/t. We have drawn 
the asymptotic configuration 0(oo, r) in Fig.2. The static equation for 0(oo, r) is 
invariant under rescaling of r. As we know from the Skyrme model, this 9(oo, r) 
is not a local minimum of energy with respect to rescaling of r since there is no 
Skyrme term in our Lagrangian. In the dynamical case under considerations, 
solutions of all different scales are acceptable for 9(oo, r). The amount of energy 
fed in by the hadron debris determines the scale of a DCC that is produced. 

From the viewpoint of total energy, the nontrivial n(x) costs nothing since 
E(x) remains at the bottom of the double-well potential and the static n(r) does 
not contribute to kinetic energy. As far as energy is concerned, the topologically 
nontrivial solutions are no different from the isospin-uniform solutions. 
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3 Isospin property 



It might appear that the topologically nontrivial DCCs have zero isospin because 
of the spherical symmetry. It is not correct. At the quantum level they are 
not eigenstates of I = but superpositions of eigenstates of many different 
isospin values. To see it, we should represent the quantum state of the classical 
configuration by the coherent state |7|, |^]: 

|n(r)0(t, r)> = exp(-if v J tr (6(t, r)n(r) • d t 7r(x) - d t 9(t, r) n(r) ■ iz(x))d 3 x^ |0), 

(12) 

where 7r(x) is the isovector quantum pion field. The exponent is invariant under 
the simultaneous isospin-orbital rotations generated by K = I + L, where L is 
orbital angular momentum. Therefore this state is an eigenstate of K = 0, not 
of / = nor I 3 = || . The exponent of Eq. (|T^) can be expressed in terms of 
the creation operators of pions with charge a = (+, — , 0), energy k(= |k|), 
and orbital angular momentum (l,m) as 

oc J + eft 1 + a^y/Uefflkdk - h.c, (13) 

where 9i(k) is the Bessel transform / \^krJiJ r i/ 2 (kr)8(0,r)r dr. Projection onto 
the N n pion state of Eq.(|l2"D is iVVth power of Eq.(|I"3"D operated on |0). By 
isospin decomposition we see that the state is not purely an isosinglet but has 
a wide distribution in /. It is remarked that the initial state of pp collision has 
K = or 1 in the case of head-on collision with zero impact parameter since 
L = and / = or 1. 



4 Momentum distribution of neutral-to-charged pion 
ratio 

We study the spectrum of the pions decaying from the topologically nontrivial 
DCCs. We focus on the correlation between the isospin and momentum distri- 
butions since it shows a distinct characteristic. Let us describe the expanding 
hadron debris of the baked Alaska scenario |10| by the isovector source p(x) 



□7r(x)). The standard method JTTJ] gives the momentum spectrum of the 
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pions with the Cartesian isospin component i and momentum k as 

where p(k) is the four-dimensional Fourier transform with k = |k| for massless 
pions, 

p(k) = J p{x)e ikx d 4 x, (15) 

and ej is a unit Cartesian isospin vector along the i-th direction. The spherically 
symmetric pion fields, 6(t,r) and n(a;) = r/r, can be produced only by p(x) of 
the same symmetry. Therefore p(x) may be expressed as 

p(x) = -p(r)5(t 2 -r 2 ). (16) 
r 

Then the k dependence of its Fourier transform is determined in the form 

p(k) = |rP(|k|). (17) 



Substituting Eq.( p!7|) in Eq.(|Tj), we find that the pion isospin is correlated with 
the momentum direction as 

2k ^ oc (k ■ e,) 2 . (18) 

Namely, neutral pions are emitted preferentially along the z-axis while charged 
pions are into the xy-p\ane. Note however that the z-axis does not necessarily 
coincide with the collision axis. For some DCC, the z-axis happens to be the 
collision axis. Then many other DCCs can exist that are related to this one by 
isospin rotations. They form one isospin family of DCC solutions. When the 
z-axis coincides with the collision axis, we expect to see a pair of parallel 7r°-rich 
bands around y\ and 2/2 in the <p-y plot, where is the azimuthal angle of the 
pion momentum k around the collision axis and y is the rapidity variable. The 
region between the two bands is filled dominantly with charged pions. When 
the isospin axis is not parallel to the collision axis, a pair of 7r°-rich domains is 
found at (0, yi) and (0 + it, y^) in the (fi-y plot. The two 7r°-rich domains are 
separated by a 7r ± -rich region (Fig. 3). 

The event-by-event distribution dP/df of the neutral pion fraction / = 
iV 7r o/(iV 7r + + N n - + N n o) can be obtained from Eq. flTif) . By integrating over the 
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angles of k, we find that the particle multiplicity is equal for all three pion charge 
states. Therefore the distribution dP/df cannot distinguish the topologically 
nontrivial DCC events from random emission events. 

We can make search of the topologically nontrivial DCCs more quantitative. 
We should first select the DCC candidates, for instance, by abundance of soft p t 
pions and find their approximate overall rest frames. To enhance the signature, 
we should select only those events with / 1/3. We then determine event 
by event the isospin polar direction z, namely the 7r°-dominant direction in the 
momentum space, by maximizing a suitably defined quantity as follows: By 
choosing a tentative z direction, compute for ir° and ^ 

C = £(k.-z) 2 , 

C± = ]T(k lX z) 2 , (19) 

where ko and k± are the unit vectors along neutral and charged pion momenta, 
respectively, and the summations are taken over all DCC pions from each event 
(N n o N n + ~ N n -). Then vary the direction of z so as to maximize the product 

C = C C ± . (20) 

The quantity C takes the maximum value when z coincides with the isospin polar 
axis. Let this maximum be C max . Then compute C for the same z, direction by 
interchanging 7r ± and 7r°. Let us call it C m i n . Then the ratio 

r — n 

b = r Tr [ ' 

^max ~r ^min 

is equal to 5/7 for the topologically nontrivial DCCs with sufficiently large N w . 
As the direction of z is varied, S sweeps between 5/7 and —5/7. In contrast, 
S is independent of the direction of z and equal to zero for the isospin-uniform 
DCCs as well as for random emission. Feasibility of this test depends on how 
large N w is. Since the statistical errors are of 0(1/ \/N^), N w k,25 will do. 
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5 Chiral rotations and disorientation 

The topologically nontrivial solutions of the form behave like ~ £ at £(= 
r/t) — > 0. At sufficiently large t, therefore, the phase pion field 9(£) approaches 
zero, that is, the background vacuum relaxes to the true vacuum at any fixed 
location off the light cone. If the DCC is defined to be the disoriented condensate 
that would approach a wrong vacuum at t — > oo, one might not want to call 
such a condensate as the DCC [[|. When we make chiral rotations on we 
obtain infinitely many more solutions that carry the same topological charge but 
different vector-isospin distributions. By actually performing the axial-isospin 
rotations, we can obtain solutions whose n(x) is not spherical at finite t. As 
t —>■ oo, these rotated solutions approach static limits with uniform isospin 
orientations at all finite locations: 

lim n(x)6(x) = n 8 , (22) 

t — >oo 

where n is the axis of an axial-vector isospin rotation (n^ = n# = n ) and 8q is 
its rotation angle {8 L = —8 R = 8 ) [fj, H]. Namely, the asymptotic disorientation 
stays nonvanishing and turns uniform. 

When 6(t,r) is not a function only of £ = r/t, we have found 8(t,r) — > r 
at r — > 0, instead of — ► r/t. It remains disoriented and nontrivial as n(x) = r/r 
even at t — oo. Chiral rotations transform n(x) = r/r into nonspherical n(x). 



6 Explicit chiral symmetry braking 

In the presence of an explicit chiral symmetry breaking, a DCC cannot really 
reach its asymptotic limit predicted by the chiral symmetric wave equations. It 
starts decaying before its kinetic energy f%8 /2 becomes comparable with the 
potential difference of the symmetry breaking m^f%(l — cos 8). However, the 
postulate of n(x) = r/r is still compatible with the explicit chiral symmetry 
breaking. Furthermore the topological charge q is defined in the same form 
and the local current conservation d^Q^ = remains valid. The only change is 
appearance of a scale breaking term proportional to m\ in the wave equation of 
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0(t,r): 

2 

+ — sin # cos # = m;; sin (23) 

With the explicitly scale dependent term present, we can no longer postulate 
= 0(0- Instead 0{x) is a function of two dimensionless variables r/m n and 
t/m^. But the wave equation still allows a spherically symmetric 0(t,r), and 
the behavior 0(t,r) ~ r at r — > is not affected by the symmetry breaking. As 
long as the form of n(r) is the same, our prediction in the <j)-y distribution of 
the decaying pions is valid with no modification. 



7 Conclusion 

We have argued that topologically nontrivial DCCs are not only possible but 
also have an equally good or bad change of being produced as the isospin- 
uniform ones. If they are actually produced, they will show a clear experimental 
signature in the <p-y plot of the decay pions. In order to produce them the 
hadron debris must also expand with a spherically symmetric isospin distribu- 
tion J oc r according to geometrical symmetry. No dynamical mechanism is 
known that suppresses such a configuration of hadron debris. In a DCC search 
the topologically nontrivial DCCs should also be searched for by the analysis 
proposed here. 
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Figure captions 



Fig.l: The scalar phase function 6(g) for a few different boundary values of 
6'(0). 6(g) oc £ is required at £ — > by the wave equation. 

Fig. 2: The asymptotic configuration of 6(t,r) at t = oo in the case that 
6(t,r) is not a function only of £ = r/i. The variable r is expressed in the unit 
of m n . 

Fig. 3: Schematic pictures of the pion charge distribution in the <p-y plot for 
a topologically nontrivial DCC whose 7 3 direction is along the collision axis (3a) 
and off the collision axis (3b). 



11 



This figure "figl-l.png" is available in "png" format from: 



http ://arXiv. org/ps/hep-ph/96073 6 1 v 1 



This figure "figl-2.png" is available in "png" format from: 



http ://arXiv. org/ps/hep-ph/96073 6 1 v 1 



